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OPTIMIZATION OF DOUBLE-CONIC INTERPLANETARY TRAJECTORIES 


by Edward A. Willis, Jr. 
Lewis Research Center 


SUMMARY 

An efficient class of high-thrust interplanetary trajectories is considered 
in which each one-way planet-to-planet transfer may consist of two distinct 
coasting arcs joined by an optimal midcourse impulse. The analysis is based on 
the successive two-body trajectory model and also assumes low circular parking 
orbits at each planet, impulsive thrust, and coplanar elliptic planet orbits. 

The location, direction, and magnitude of the midcourse impulse are optimized 
directly by a convenient numerical search procedure. 

These double-conic trajectories are compared with conventional single-conic 
trajectories using one-way and stopover round trips to Mars as typical examples 
with total velocity .increment, trip time, and stay time for criteria. It is 
shown that double-conic trajectories significantly reduce the total velocity in- 
crement for one-way transfers which cover heliocentric travel angles greater than 
4 radians. Round-trip durations of 300 to 1000 days with stay times of 0 to 
140 and 450 days are then studied in 1971 and 1980. A 500-day trip with a 40- 
day stay time is also investigated over a range of synodic periods from 1971 to 
1980. The double-conic trajectories yield attractive reductions in the propul- 
sive velocity increment sum, and/or lower Earth-approach (reentry) veloc- 
ities, for round-trip times of 360 to 700 days. These advantages are available 
in all synodic periods, but are most pronounced in the more difficult periods. 


INTRODUCTION 

Trajectory studies are an essential part of interplanetary mission analysis. 
Trajectory data are not only used directly in the mission calculations, but cri- 
teria such as the total velocity increment ^ AV, total trip time T^ 0 ^., and 

stay time T s by themselves provide valuable guidelines for selecting candidate 
trajectories. For instance, trajectories with minimum ^ AV often yield low 
if not minimum space vehicle weight. 

Previous surveys of high-thrust trajectories, such as references 1 to 3, 
have assumed that major propulsive maneuvers only occur near a planet terminal. 
Each one-way heliocentric transfer would then consist of a single coasting arc, 
and a round trip would consist of two such single-arc transfers. This assump- 
tion leads to a convenient numerical computational procedure but not always to 
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minimum total velocity increments. Rigorous investigations, such as reference 
4, have shown that optimal one-way transfers in some cases consist of two 
coasting arcs joined by a midcourse impulse rather than one single arc. This 
raises the possibility, not considered in reference 4, that double-arc trans- 
fers might also be useful as one or both legs of a round trip. The rigorous 
variational solutions are valuable as standards of comparison but have not re- 
ceived a wide practical application to high-thrust trajectories because of in- 
herent computational difficulties. 

The two objects of this report, therefore, are to develop a more effective 

technique to find optimal double-arc transfers, and to present examples illus- 

trating the characteristics and utility of this class of trajectories for both 
one-way transfers and round trips. 

The following analysis considers trajectories in which each one-way planet- 
to-planet transfer may consist of two distinct coasting arcs joined by a mid- 
course impulse. These are referred to as "double-conic" or "three-impulse" 
transfers in contrast to the conventional single-conic transfers assumed in ref- 
erences 1 to 3. The two arcs lie in a single plane as opposed to the "broken- 

plane" transfers discussed in reference 5. The magnitude, direction, and loca- 
tion of the midcourse impulse are optimized with the aid of a direct method 
(first described in ref. 6) rather than by a variational technique. The present 
formulation includes single-conic transfers as a special case in which the mid- 
course impulse vanishes. Also included as a special case is the "perihelion- 
propulsion" trajectory in which the midcourse impulse occurs at the common per- 
ihelion of the two arcs. All previously established optimal results such as the 
classical Hohmann transfer and the 360° three-impulse transfers discussed in 
reference 4 are recovered when appropriate boundary conditions are applied. The 
resulting program is not only computationally efficient, but is also very flex- 
ible in that various criteria of merit, boundary conditions, and constraints can 
be easily accommodated. 

A family of potentially desirable Earth-Mars trips is studied as a typical 
application of this analysis; the velocity increment sum is used as an illustra- 
tive criterion of merit. Round trips lasting 300 to 1000 days, with stay times 
at Mars of 0 to 140 and 450 days, are considered in the oppositions of 1971 and 
1980. A 500-day round trip with a 40-day stay is evaluated in all oppositions 
from 1971 to 1986. The effects of atmospheric braking and trajectory profile 
are also discussed, and double-conic trajectories are compared with conventional 
and Venus -swingby trajectories. 

SYMBOLS 

C coefficient defined in appendix A 
D date, Julian day 
e eccentricity 
M planetary motion, rad/day 
N revolution index 
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T time, days 

AT terminal -to -terminal travel time, days 

5T time error, defined in appendix A 

V velocity, miles/sec 

AV n impulsive velocity increment at terminal n, miles/sec 



total velocity increment for n impulses, miles/sec 


(X) 

“P 


denotes time-average of X 

planet perihelion position (see table V) 


p heliocentric travel angle (see fig. 4), rad 

T] trajectory true anomaly at planet terminal, measured positive in 

direction from midcourse maneuver toward planet (see fig. 4, views 
A and B) , rad 

0 orbital true anomaly of planet, measured positive, counterclockwise, 

rad 

K travel angle distribution parameter, (ij/p 

(j. partial travel angle (see fig. 4), rad 

v trajectory true anomaly at midcourse maneuver, measured positive to- 

ward planet (see fig. 4, views A and B), rad 

p heliocentric radius, astronomical units, or planetocentric radius in 

planet radii 

T planet orbital period, days 

f trajectory path angle (fig. 4), rad 

Subscripts : 

e atmospheric entry 

1 inner 

inw inward 

k kick point (i.e., location of midcourse maneuver) 
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o outer 

out outward 

opp opposition 

p perihelion 

pi planet 

po parking orbit 

r relative 

s stay 

si sphere of influence 

t trajectory 

tot total 

00 denotes conditions at great distance from center of gravitational force 

© Earth 

cf Mars 

1 Earth departure 

2 outbound kick point 

3 Mars arrival 

4 Mars departure 

5 inbound kick point 

6 Earth arrival 

ANALYSIS 

Several major steps are involved in seeking desirable trajectories: 

(1) Criteria of merit must be defined. 

(2) An appropriate computational model must be devised. 

(3) A suitable optimization technique must be applied to this model to ob 
tain the best possible values of the selected criteria. 

These steps will be developed below for one-way and round-trip trajectories 
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using both single- and double-conic transfers. 


Criteria of Merit 

Many independent variables must be considered in evaluating final mission 
criteria such as space vehicle weight or cost. Several useful mission param- 
eters, however, may be studied on the basis of trajectory information alone. 

These include the following: 

(1) The total velocity increment ^ AV 

(2) The total trip time T^ ^ 

(3) For round trips, the stay time T g at the destination planet 

The most important of these is the total velocity increment. The mass 
ratio required for each propulsive maneuver is directly related to the. associ- 
ated velocity increment, and the sum of these increments is certainly the best 
single indication of space vehicle weight that can be obtained without making 
numerous assumptions about the vehicle systems, payloads, and mission objectives. 
The velocity increment sum is therefore used as the primary criterion of merit 
for the purpose of developing and illustrating the trajectory optimization pro- 
cedure. It could easily be replaced by cost, weight, or other criteria in 
order to study a specific mission or program; this substitution would not in- 
volve any essential change in the method of analysis. 

The trip and stay times are treated as secondary criteria. A low value of 
T^ 0 t would indicate a minimal exposure to the hazards of space, and for round 
trips, the stay time limits the time available to accomplish the primary mission 
objectives . 

With these criteria, the optimization procedure to be developed will yield 
minimum ^ AV one-way transfers which satisfy specified boundary conditions 
such as travel time and heliocentric travel angle. These optimal one-way trans- 
fers may be of interest for probe missions, and are also used to construct min- 
imum ^ AV round trips with prescribed trip time T^^ and stay time T s . 


Appr oximat ions 

This analysis is based on the impulsive two-body trajectory model as de- 
scribed in reference 1; that is, propulsive maneuvers are treated as impulses 
and the actual n-body problem is replaced by a sequence of two-body coasting 
arcs. Successive heliocentric arcs are joined by a midcourse impulse; helio- 
centric and planetocentric arcs are related by matching conditions applied at 
the "sphere of influence." Inside the sphere of influence, only the planet's 
gravitational field is considered, and outside, only the Sun's. The spheres of 
influence are taken to be of negligible size compared to interplanetary dis- 
tances, but much larger than the parking orbit radii. It is assumed that trans- 
fers begin and end in circular planetocentric parking orbits, and that the plan- 
ets themselves lie in elliptic^ coplanar heliocentric orbits. As shown in ref- 

TZ 

Circular planet orbits are temporarily used for illustrative purposes in 
this section, but the data to be presented in the RESULTS AND DISCUSSION are 
based on elliptic orbits. 
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erence 5, the ^ AV data for optimal three-dimensional (broken -plane) trans- 
fers are well approximated by coplanar elliptic data; this will be discussed 
later (p. 19) . 


Under these assumptions, the well-known conic orbit equation (derived in 
ref. 7, for example) is applicable. The form appropriate to the present anal- 
ysis is 


Pp(l + e ) 

1 + e cos 0 


(i) 


Computation of Single-Conic Transfers 

Single-conic transfers merit individual attention because they occur fre- 
quently as a special case of the more complex double-conic transfers and have 
several elements in common with them. 

The heliocentric portion of a typical one-way single-conic transfer is il- 
lustrated in figure l(a). The space vehicle departs from the inner planet by 
means of an impulse applied at the circular parking orbit. It pierces the 
sphere of influence at point 1 and follows the heliocentric arc 1-3, describing 
a central travel angle P ou t> and - finally encounters the destination planet's 
sphere of influence (point 3) after the travel time ATq.j has elapsed. The 
transition from heliocentric to planetocentric coordinates is illustrated in 
figure l(b). The vehicle follows the planet -approach hyperbola 3-3' from the 
sphere of influence to the perigee, which is at the desired parking orbit ra- 
dius. The impulsive thrust AV changes the vehicle's velocity from Vp to 
the local satellite velocity V po . 


Heliocentric trajectory . - The planet positions and hence the trajectory 
terminal conditions p, p 0 , and pj_ corresponding to a given travel time and 
planet encounter date may be computed from elliptic ephemeris data as described 
in appendix A. These quantities together with the travel time AT^_ 0 com- 
pletely specify the elements of the transfer arc. 


The necessary relations are derivable from the basic equation (l). When a 
trial value r] i is used for the trajectory true anomaly at the inner terminal, 
the eccentricity is 


_ Pq ~ Pj 

e t " Pi COS Tl i - P Q COS T) 0 

where the true anomaly at the outer terminal is 

Tj = T). + P 
'o *1 


( 2 ) 


(3) 
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The trajectory perihelion radius is then 


1 + e + cos ru* 

= <* "U 


( 4 ) 


With equations (2) and (4), the required value of may be determined by 

iterative solution of the time equation (illustrated here for elliptic arcs) : 



(The iteration cycle is described in appendix A.) The velocity and path angles 
for both the trajectory and the planet orbits may then be computed at each en- 
counter from the expressions 


tan i|r = 


e sin 9 
1 + e cos 9 


V = 18.5058 


2 

P 



(Equations (5) to (7) are also derived in ref. 7.) 


( 6 ) 

(V) 


Planetocentric trajectory . - The transition from heliocentric to planeto- 
centric coordinates is illustrated in figure l(b) . Because of the assumption 
made about the size of the sphere of influence, the planetocentric hyperbolic 
excess velocity may be computed by evaluating the law of cosines at the 

point of intersection of the trajectory and the planet orbit: 


V 2 = V? + V 2 - 2V, V . cos (ty , - If ,) ( ? ) 

CO t pi t pi t pi 

By the two-body assumption, is a constant of the planetocentric motion and 

may be used in the energy equation to evaluate the impulsive velocity increment 
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Mars orbit 




Figure 1. - Typical interplanetary transfer. 



Travel angle, ft rad 

Figure 2. - Characteristics of single-conic Earth-Mars 
transfers. 


required to enter a circular parking 
orbit whose circular velocity is V : 



(Values of V corresponding to 
p^ o = 1.1 are listed in table V, ap- 
pendix A. ) 

Optimization of Single-Conic 
Trajectories 

One-way transfers . - For a given 
travel time, the criterion function to 
be minimized is 

VAV = AV . . + AV _ (10) 

pl,i pl,o 


It is important to note, from the foregoing analysis, that a transfer of this 
kind can be optimized only with respect to two parameters, the travel time AT 
and the travel angle p. (Note that in a given synodic period, and with given 
AT, the choice of the travel angle (3 is equivalent to choosing the planet en- 
counter date.) The influence of these parameters on ^2 AV is illustrated in 
figure 2 for Earth -Mars transfers in which circular -coplanar planet orbits were 
assumed for simplicity. The ^2 AV is plotted against p for AT^^ = 250, 
350, and 450 days. If both P and AT 0 _^, are left open, optimization yields 
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the well-known Hohmann type of trans- 
fer. These have a low ^ AV but in- 
volve definite values of 3 and AT, 
which cannot always be incorporated in- 
to a desirable round trip. If one pa- 
rameter (e.g., AT) is fixed in advance, 
then optimization can be only carried 
out with respect to the remaining one 
(e.g., 3). If both parameters are 
fixed, then the transfer is completely 
specified and not subject to any kind 
of optimization. Note in particular 
the very high AV T s that occur for 
large values of 3 . 

Round trips . - A round trip con- 
sists of two one-way transfers separated 
by a stopover at the destination planet, 
as illustrated in figure 3. The com- 
plete optimization of a round trip using 
single-conic transfers would therefore 
require the determination of the most 
advantageous travel time and angle for each leg. These must be determined in 
consistency with the Earth rendezvous condition. That is, since round trips 
begin and end at Earth, the total angular motion of the space craft during the 
total trip time mus ^ exactly match (within a whole revolution) the 

Earth* s motion during the same time interval: 

Pout + Pinv - <V T tot - <V T S + 2 ” K < N - °- ±2 > • • ■ ) 

AT out + &T im, * T tot - T s < 12 > 


Figure 3. - Typical high-thrust trajectory for Mars round trip. Total 
trip time, 300 days; stay time, 40 days. 


As the rendezvous condition only needs to be satisfied within whole revolutions, 
the revolution index N is also subject to choice. A distinct family of round 
trips corresponds to each choice of N; values of 0 and -1 will be of in- 
terest for the Mars trip discussed later (see p. 23.) 

"Double -Hohmann" round trips in which the travel times and angles of the 
outbound and return transfers are optimized independently yield low ^ AV 
but require specific values of T.^^ an( ^- When T^ o ^ and T g take on 

arbitrary values, only two parameters remain free for optimization. Whatever 
two are used (for instance, the planet encounter date and the outward travel 
time), the problem ultimately reduces to choosing the optimum distributions of 
the available travel angle and travel time (cf. eqs. (ll) and ( 12 )) between 
the outward and return legs . 

As a highly simplified example, consider a 500-day Earth-Mars round trip 
with negligible stay time and N = 0. (Circular-coplanar planet orbits are 
still assumed; fig. 2 gives the ^ AV as a function of 3 and AT for both 
outward and return transfers.) Then, according to equations (ll) and (12), the 
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two travel angles P Qut and P inw must sum to about 500° or 8.7 radians; the 

two travel times ^- ou t an< ^ ^inw mus ^ U P ' t:o 500 days. A natural first 

trial in seeking the optimum distribution for this trip would be to use the 
Hohmann outward transfer; this requires ^out ^ .14 radians , AT^^. = 250 days, 

and ^ AV out = 3.4 miles per second as shown on figure 2. The return transfer 
must then have AT inw *» 250 days and 3j_ nv = 5.56 radians. With such a long 
travel angle, a ^ ^inw = ^.l miles per second is required for the return 

transfer. The grand total for the round trip is then 16.5 miles per second. By 
similarly evaluating other distributions of travel angle and time, the best 
round trip is found to consist of symmetric outward and return legs, each having 
a travel time of 250 days, a travel angle of 4.35 radians, and a £ AV of 
6.8 miles per second. The grand total AV for this round trip is then 13.6 
miles per second. Neither the outbound nor the return transfer in this example 
has an individually optimum combination of 3 and AT, yet together they com- 
prise an optimum round trip. 

To summarize the foregoing remarks, a single-conic one-way transfer may be 
optimized only by varying its boundary conditions, that is, 3 and AT. If the 
transfer is supposed to be part of a round trip, even the possible choices of 
3 and AT are restricted by the Earth rendezvous condition. 


Computation of Double-Conic Transfers 

Double-conic transfers are more responsive to optimization than are the 
single-conics considered previously. Not only the terminal conditions (3 and 
AT) but also the parameters describing the location, magnitude, and direction 
of the midcourse impulse can be varied. Because of this added flexibility, a 
double-conic transfer can be optimized to a significant extent even when 3 
and AT are fixed. 

This section will develop the procedure for computing the elements of a 
double-conic transfer for arbitrary terminal conditions and an arbitrary mid- 
course impulse. The following section will then consider optimization of double- 
conic transfers as applied to both one-way and round trips. 

Geometry of double-conic transfers . - Figure 4 illustrates the geometry of 
a one-way double-conic transfer. As was the case for single-conic transfers, 
the terminal radii p Q and p. and the travel angle 3 between the terminals 

are fixed when the planet encounter date D Q and the travel time AT are 
given. The kick impulse AV^ occurs at radius p^ and divides the transfer 

into two conic arcs. The arc from the inner terminal covers a central angle 
P-p and it has true anomaly Vj_ and path angle ^ at the kick point. Like- 
wise, the outer arc describes a central angle of p 0 = 3 - and has true 
- anomaly and path angle at the kick point. 

o K, o 

Heliocentric trajectory . - The elements of each of the two conic arcs which 
comprise the heliocentric transfer are uniquely determined by the geometrical 
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View A: Inner trajectory 
segment. 


Figure 4. - Geometry of Mars-Earth double-conic trajectory. 


quantities p . , u . , p . (i . p, , and , which were illustrated previ- 

ously. The necessary relations may be easily derived from the basic equations 
(l) and (6). The trajectory true anomaly at the inner terminal is 

\ + v ± (13) 


The eccentricity of the inner arc is then found in terms of v-j_ by evaluating 
equation (l) at the kick point and at the inner terminal and by eliminating 

between the resulting two expressions: 


e 


i 


Pi - Pk 

p k cos v ± - p i cos rj i 


(14) 


Substituting equation (l4) into equation (6), evaluating the result at the kick 
point, and solving for yield 


v . 
1 


tan 


-1 


tan t k ^1 - cos q i ) 


1 - 


Pi 


sin tan \|r k ^ 


rad 


(15) 


when the identity for cos(x + y) is used. Finally, the perihelion radius of 
the inner arc is 
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P k (! + cos v ± ) 
p p,i 1 + e j _ 

The required procedure is then to 
first compute v^, then e i , and 

finally p . . The time of transit 

along the inner arc may then he 
computed directly from equation (5) 
evaluated between and t)^. 

These results apply also to 
the outer subarc (replace subscript 
i by o) and have the same form 
whether the original transfer is 
inward or outward. This follows 
from the sign convention, illus- 
trated in figure 4, that B and 
3 5 7 9 11 13 15 ^ * 

Kick radius, p k . astronomical units ' ' H are always positive and that 

true anomalies range from -it to 

Figure 5. - Mars-Earth transit times for double-conic transfers. , , . , . , , 

n and are measured positive toward 

the planet associated with each 

subarc, that is, away from the kick point. (The sign convention for the path an- 
gles an<i o is in(3 - uce( i by means of eq. (6).) The total travel time 

for the entire transfer may then be computed by applying equation (5) to both 
subarcs : 

’ll n 0 

+ AT o (17) 

V i 

Equations (13) to (17 ) specify the transfer in terms of the four independ- 
ent parameters p^, p^., ^ and o’ Th ese j together with the terminal 
conditions, also implicitly define the midcourse impulse. 

One independent parameter must be used to satisfy the travel time con- 
straint. An advantageous choice for this purpose is the kick radius p^. With 
given values of 3, p^, p^, p Q , and ^ the travel time (eq. (17)) 

varies monotonically and almost linearly with p k as illustrated in figure 5. 
Although pk cannot be computed in closed form in terms of AT^_ q , the required 

value may be determined from a simple iteration process of the type described in 
appendix A. For Earth -Mars transfers, a good initial guess for p^ may be ob- 
tained from figure 5. 

The preceding derivation gives the elements e, pp, and v of a double- 
conic trajectory which satisfies prescribed terminal conditions and travel time 
constraint in terms of three independent parameters p^, ^ ±> and tk,o* Th ese 

parameters implicitly define the midcourse impulse. The magnitude of the mid- 
course impulse may then be obtained from these elements by evaluating equations 
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(6) and (7) at the kick point for both the inner and outer arcs and substituting 
the resulting velocities and path angles into the law of cosines: 

AV v = + [v 2 + V 2 - 2V.V cos (\|r, . + t. )1 ^ (18) 

k i o l o ' T k,i k,o'l ' 

Finally, the planetocentric velocity increments are computed from equation (9), 
just as they were for single-conic transfers. 


Optimization of Double-Conic Trajectories 

For one-way transfers, the criterion to be minimized is 

L AV = AV . . + AV + AV (19) 

pl,i k pi, o 


For round trips the object is to minimize the grand total AV; that is, 


Z AV + + - Z AV * + E av. 

tot out mw 


( 20 ) 


One-way transfer, fixed terminals . - The velocity increment sum for a 
double-conic transfer can be minimized even when the terminal conditions (de- 
fined by p , p^, 3, and AT) are prescribed in advance by making use of the 
extra degrees of freedom associated with the midcourse impulse. The problem is 
simply to choose the three parameters q^, V k i , and V k Q which control the 

midcourse impulse so that equation (l9) is minimized. An equivalent but more 
convenient choice of control variables is obtained by the transformation 


A 


Pi = 0K 

° < K < 1 

*k,i “ *k 

- 1 < < 1 l 

2 “ k ~ 2 

= 'if 

k,o r 

*k'*k"£ <*r<*k + ^ 


( 21 ) 


Here K (= q^/p) is just the angular location of the kick expressed as a frac- 
tion of 3, T|r k is the path angle of the inner arc at the kick point, and 

is the discontinuous change in the path angle in going from the inner arc to 
the outer. The restrictions on the ranges of the control variables are needed 
to exclude retrograde arcs . 


It should be noted that the control variables are bounded. Thus, an ex- 
haustive investigation of the control space may be made in the search for mini- 
mum ^ AV. This is an important point since the existence of multiple local 
minima cannot be ruled out by a priori considerations. A typical example of an 
exhaustive search is shown in figure 6. Contours of constant ^ AV are 
plotted on the plane = 0 of the control space (k, i k , i r ) for a 5.5-radian, 
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250-day, Mars -Earth transfer. The plane 
= 0 is the locus of double-conic 
transfers using a tangential midcourse im- 
pulse. It turns out that these "smooth 11 
transfers are, in fact, optimal for Earth- 
Mars transfers of practical interest. The 
heavy curve (^jAV = 13 miles/sec) forming 
the left and lower boundary is the trace 
of the single-conic transfer (i.e., where 
AV^. = 0) . This curve is simply a plot of 

path angle against dimensionless central 
angle for a heliocentric elliptic arc 
that traverses 5.5 radians in 250 days. 

The curve forming the upper and right- 
hand boundary represents a family of 
double-conic trajectories that also re- 
quire 13.0 miles per second. These two 
curves form a closed contour bounding the 
region of control space in which double- 
conic trajectories are more efficient 
than single conics. The unique minimum 
of 53 AV, 8.9 miles per second, occurs 
near the center of this contour system where = 0.033 radian and K = 0.454. 

These values are numerically close to those of the single-conic perihelion for 
which ^ = 0.0 and K = 0.41, but they decrease the 53 AV by 4.1 miles per 

second or 30 percent. 

Search procedure . - Since the contours of constant 53 AV cannot be rep- 
resented mathematically in closed form, an efficient numerical search procedure 
is needed to find the optimal kick impulse controls for the many transfers of 
interest. One suitable procedure is developed subsequently using the two- 
dimensional case where = 0 as an example. (The same basic procedure can 

also be used for a three-dimensional search where K , and \|/ are all var- 

ied.) 
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Figure 6. - Velocity-increment contours for 5.5-radian, 
250-day, Earth -Mars transfer. 


In the neighborhood of the minimum 53 AV point, the contour system forms 
a long, narrow, and fairly straight valley. An efficient numerical search pro- 
cedure for such contours is the method of parallel tangents (PARTAN) as devel- 
oped in reference 8. A unidimensional search along an arbitrary line, such as 
a - b in figure 6, yields a relative minimum at point e. A second search along 
a line c - d (parallel to a - b) yields another relative minimum at f. A third 
search along the line through e and f yields a candidate minimum at g. If 
the contours in the neighborhood of the minimum are concentric ellipses, then, 
as proved in reference 8, point g is the exact center of the contour system, 
and J3 AVg is the true minimum. For more general contours, the actual center 
can be approached within arbitrary tolerances by iterating this procedure 
starting with a new arbitrary line through g. 

Starting the search. - For transfers where p < 2it radians, the single- 
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conic trajectory may be used to furnish starting conditions that ensure rapid 
convergence of the PARTAE search cycle. It may be observed from figure 6 that 
the perihelion, tyj-, = 0, is "closer" to the minimum than any other point along 
the single-conic boundary; this is a natural starting point for the search. The 
elements of the single-conic transfer may of course be determined as in the pre- 
ceding section. The coordinates of its perihelion, in the control space, are 


K 


P 


3 



( 22 ) 


V * 0 
r ' P 

This choice is not possible if the single-conic transfer does not pass through 
its perihelion. Reasonable initial values then are K = 0.5, Tjr = 0, and ^ 

from equation (6) . These choices together ensure rapid convergence and guaran- 
tee that the search will yield single-conic transfers when those are optimal. 


Although not strictly necessary, it is advantageous to take the initial 
line a - b in the steep descent direction, that is, in the direction of -grad 
( ^ AV) . Then the major part of the saving occurs early in the PARTAE cycle - 
a time-saving feature for preliminary surveys. 


This gradient is perpendicular to the single-conic contour; hence, the 
slope of a - b for steep descent is 


ak 


-1 


b.-b 


dj 

d/C 


(23) 


p, single conic 


The necessary derivative may be obtained directly from equation (6) upon noting 
that 


_d 

d/C 


TZ * 3 


dp j _ 



(24) 


That is, 


d/C 


-e(3 cos 


cos Vn + e 


(l + e cos v . )' 


(25) 


When equation (25) is evaluated at the perihelion > v/here = 0, equation 

(23) becomes 


A \ 



1 + e 

ep 


(26) 
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This is the steep-descent slope for the initial line a - b, where point a (the ' 
single-conic perihelion) is defined by equation (22) . 

It is considerably more difficult to obtain satisfactory convergence for 
transfers where (3 > 2« radians, since a single-conic transfer cannot be used 
to obtain a good initial starting point and direction. Such transfers, however, 
appear to be of little more than academic interest, at least for Mars trips, 
and will not be further discussed. 


Convergence . - At the end of each PARTAN cycle, the candidate minimum at 
g must be further investigated to be sure that it is actually a minimum and not 
a saddle or inflection point. A necessary and sufficient condition that y AV 
be minimum at g is that 


df 

5c 


df 

“ ^ 

d 2 f 

d/c 2 


= 0 


> 0 


and 


a 2 f a 2 f 

dK 2 d ^ 2 


where 


£av = f ( K ,* k ) = f g + & 


AK + 


df 


A *k 



2 d 2 f 
^ + 2 


*\ 



> 0 


(27) 


J 


-v 2 

S f 

+ ry 


Af, 


(28) 


g 


The necessary partials may be numerically approximated, by the usual algebraic 
means, after computing five test points in a small neighborhood of g. If the 
conditions of equation (27) are satisfied within acceptable tolerances, the 
search is terminated at g; if not, it is iterated until they are by using a 


new initial line through g. 

It is convenient 

in the steep descent direction also, since the 


Sf 

A *k 


A/C 

as 

df 


(29) 


initial, 
new line 


cty 
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Figure 7. - Comparison of single- and double-conic Earth-Mars transfers, 
and the required partial derivatives are already available in equation (27). 

One-way transfers; variable terminals . - The preceding discussion consid- 
ered the optimization of one-way transfers with fixed terminals defined by 3 
and AT. Further reductions in ^ AV may be expected for double-conic 
transfers in those cases where the terminal conditions (3 and AT) may be op- 
timized in addition to the midcourse impulse. As was the case for single-conic 
transfers, an optimum value of (3 occurs for each choice of AT; and if both 
3 and AT are left open, an overall minimum ^ AV exists for one specific 
combination of these parameters. These facts, plus a general comparison of 
single- and double-conic transfers to Mars, are illustrated in figure 7. Again, 
circular coplanar planet orbits were assumed for simplicity with = 1.53 as- 
tronomical units; the AV sum is plotted against the travel angle (3 for 
travel times of 250, 350, and 450 days. The single-conic transfers indicated by 
dashed lines are the same ones that were shown before in figure 2. 

Double-conic transfers, denoted by the solid curves, show a sizable advan- 
tage compared to single conics in the long travel-angle region. The reduction 
in ^2 AV becomes noticeable at (3 « 4 radians, reaches a maximum for 3 * 2 tc 
radians, and extends out to p * 7 or 8 radians. If 0 is outside these lim- 
its, the results of the double-conic optimization procedure reduce to single- 
conic transfers. For instance, the optimum AV^. approaches 0 continuously as 

3 approaches about 4 radians from above. 

For long travel times (over 350 days) a new minimum value of ^ AV, lower 
than the single-conic minimum, occurs in the long-angle region. The correspond- 
ing optimum value of 0 is also greater than it was for the single-conic trans- 
fer. For shorter travel times (under 350 days), the optimum trajectories reduce 
to single-conics before the best value of 0 is reached. In particular, when 
both 0 and AT are left open, the classical single-conic Hohmann transfer is 
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recovered. 


In summary, the two characteristic properties of double-conic transfers are 
as follows: (] T ) For all travel times, the optimal double-conic transfers re- 
quire significantly lower AV than the corresponding single-conic transfer, 

if the travel angle 0 is required to be larger than about 4 radians, and (2) 
for shorter angles, the kick impulse AV k vanishes, and the optimal double- 
conic transfers reduce to single conics. 

From these observations, it may be anticipated that double-conic transfers 
will be of the greatest utility where external constraints, such as the Earth 
rendezvous condition for round trips, preclude the selection of both optimum 0 
and AT. 


Round trips * - A round trip in which a double conic is used for either the 
inward or outward transfer, or both, is referred to herein as a double-conic 
round trip. Two main steps are required to optimize a round trip in which 
double-conic transfers may be used: (l) For each trial distribution of 0 and 

AT, both one-way legs must be optimized as described previously, and (2) the 
optimal distribution of 0 and AT must then be found in consistency with the 
Earth rendezvous condition. As was the case when single-conic transfers were 
used, the rendezvous condition (eqs. (ll) and (l2)) implies that a round trip 
with arbitrary values of and T s , and N = 0, must include at least one 

long-angle transfer. It is of interest to reconsider the 500-day, zero-stay- 
time Mars round trip that was used as an example previously. Again, the two 
travel angles must total 8.7 radians; beginning with the Hohmann outward trans- 
fer (0 = 3.14 rad, AT QUt = 250 days, and £ AV out = 3 - 4 miles /sec), the re- 
turn transfer must have ATj_ nw = 250 days and 0 inw =5.56 radians. From fig- 
ure 7, the ^ AV"in V is ^en miles per second for a grand total of 12.3 

miles per second. This is already a substantial reduction (l.3 miles/sec) over 
the minimal single-conic round trip for which ^ A ^tot was 43 m i 4es P er 

second. It is of interest to note that if double-conic transfers had been ap- 
plied with the same time and angle distribution that were found optimal for the 
single-conic trip (0 = 4.35 rad and AT = 250 days for each leg), only a 
small improvement would have been found. This simple and incomplete example is 
sufficient to illustrate (l) that reoptimization of the travel times and angles 
is required to develop the full potential advantage of double-conic trajector- 
ies, and (2) that the resulting optimal distributions are considerably differ- 
ent than those found for single-conic trajectories. 


Computational Technique 

The procedures and equations indicated above were programmed for numerical 
solution on a high-speed digital computer. Planet orbit-element and physical 
data were taken from reference 3. The resulting program may be used for the 
computation and optimization of single- and/or double-conic one way and round 
trips between any combination of planets. Optimization criteria other than 

AV, such as space vehicle weight or cost, could be accommodated without re- 
quiring a basic reformulation of the method. 
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RESULTS ABED DISCUSSION 


Families of one-way transfers and round trips to Mars were studied as a 
typical application of the present analysis. The illustrative examples dis- 
cussed in the preceding section assumed circular planet orbits, but planet-orbit 
eccentricities are accounted for in the data presented here. The results illus- 
trate the main features of double-conic trajectories and may be directly com- 
pared with previous results (refs. 1 to 3), which were based on single-conic 
transfers . 


One-Way Transfers 

One-way transfers are of direct interest for probe missions and are also 
the basis for constructing round trips. 

A spectrum of potentially desirable outward and return transfers in 1979 to 
1980 is presented in tables I and II. The 1980 synodic period was chosen for 
discussion because it is one of the more "difficult" periods and is among the 
earliest in which a manned round-trip mission to Mars might be contemplated. 

The data shown were computed on the basis of minimum ^ AV and will be 
used later to construct round trips with minimum ^ AV. These data can also 
be used for preliminary gross weight and other mission calculations. 

Table I lists the parameters of interst for single-conic outward transfers, 
that is, ^ AV, the individual velocity increments, the travel angle, and the 
minimum solar approach radius, for travel times ranging from 140 to 280 days, 
and for Mars encounter dates ranging from Julian day 2444380 to 2444460. The 
corresponding parameters are listed in table II for double-conic inward trans- 
fers with travel times of 220 to 270 days and Mars encounter dates from 2444420 
to 2444500. The data in table I can also, with little error, be associated with 
certain inward transfers, and the table II data associated with a group of out- 
ward transfers. The reasons why this is true and the procedure for doing so are 
discussed in appendix B. Direct use of the tables I and II data will result in 
a family of "short-long" round trips in which the long-angle double-conic leg 
occurs last. The "long-short" round trips associated with the alternate inter- 
pretation then use the long-angle double-conic transfers as the outbound leg. 
These two types of round trips are compared later (see p. 29) • 

Effect of planet-orbit inclination . - The outward transfers listed in 
table I involve short travel angles, for which single-conic transfers are op- 
timal. These data may therefore be compared directly with the results of ref- 
erences 2 and 5, in which inclination effects are considered. This comparison 
is presented in figure 8, in which AV is plotted against Mars encounter 

date (and travel angle) for an Earth-to-Mars travel time of 200 days. The 
coplanar -elliptic results of the present analysis are shown by the solid curve. 
The dashed curve indicates the results of reference 2, in which a three- 
dimensional trajectory is assumed to lie in a single transfer plane. These 
correspond very closely to the present results except near a travel angle of 
n radians. In this region, the single transfer plane must be inclined nearly 
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CHARACTERISTICS OP OPTIMAL EARTH-TO-MARS TRAJECTORIES FOR 1980 


Total velocity Earth departure Kick, Mars arrival 

increment, a V 

V w1 , 2* AV 3 , V , 

miles/sec miles/sec miles/sec ^Hes/sec mlle3 / 3ec mjl " a / Se 


Minimum 

radius, 

p min' 

astronomical 

units 


Central Transit 
angle, (3, time, T e-C ?, 
rad days 


Julian date 
at Mars 
arrival 


6.861 

3.190 

4.256 

__ 

3.667 

4.927 

0 . 983 

1.663 

140 

6.427 

2.813 

3.510 

-- 

3.611 

4.802 

.984 

1.841 

150 

6.136 

2.578 

2 . 977 

— 

3.555 

4.797 

.985 

2.018 

160 

5.945 

2.454 

2.660 

-- 

3.488 

4.719 

. 986 

2.195 

170 

5.822 

2.412 

2.545 

-- 

3.408 

4.625 

.988 

2.372 

180 

5.738 

2.428 

2.589 

-- 

3. 308 

4.055 

.986 

2.548 

190 

5.675 

2.484 

2.738 

-- 

3.189 ! 

^.363 

.980 

2.722 

200 

5.623 

2.506 

2.948 

-- 

3.054 1 

4.201 

.972 

2.897 

210 

5.578 

2.669 

3.191 


2.907 

4.071 

. 962 

3.070 

220 

6.901 

3.686 

5.116 

-- 

3.212 

4.392 

0.983 

1.565 

140 

6.313 

3.154 

4.190 

-- 

3.156 

4.325 

. 983 

1 . 742 

150 

5.906 

2.783 

3.445 

-- 

3.120 

4.281 

.984 

1.920 

160 

5.633 

2.543 

2.890 

-- 

3.087 

4.241 

. 985 

2.098 

170 

5.434 

2.409 

2.536 

-- 

3.042 

4.186 

. 986 

2.275 

180 

5.342 

2.354 

2.378 

-- 

2.985 

4.117 

. 988 

2.451 

190 

5.271 

2.358 

2.389 

-- 

2.911 

4.025 

.987 

2.627 

200 

5.223 

2.403 

2.520 

-- 

2.816 

3.908 

. 983 

2.802 

210 

5.188 

2.478 

2.724 


2 . 707 

3.771 

. 976 

2.976 

220 

6.426 

3.644 

5.048 i 


2.779 

3.862 

' 0.983 

1.645 

150 

5.869 

3.132 

4.149 ! 


2.734 

3.805 

1 . 983 

1.823 

160 

5.480 

2.766 

3.410 


2.711 

3.777 

.984 

2.001 

170 

5.218 

2.524 

2.841 

-- 

2.6 92 

j 3.753 

. 985 

2.178 

180 

5.048 

2.380 

2.455 

-- 

2.665 

i 3.718 

. 986 

2.356 

190 

4.940 

2.314 

2.258 

-- 

2.623 

i 3.665 

. 988 

2.532 

200 

4.878 

2.307 

2.235 

-- 

2.508 

i 3.594 

. 98 9 

2.708 

■ 210 

4.840 

2.343 

2.345 

-- 

2.495 

3.500 

.985 

2.683 

220 

4.822 

2.411 

2.541 


2.407 

3.386 

. 978 

3.057 

230 

6.041 

3.607 

4.987 

__ 

2.431 

3 . 417 

0 . 983 

1.727 

160 

5.508 

3.114 

4.115 

-- 

2.391 

3.364 

. 983 

1.905 

170 

5.133 

2.756 

3. 387 

-- 

2.374 

3.342 

.984 

2.083 

180 

4.872 

2.512 

2.811 

-- 

2.364 

3.328 

. 985 

2.261 

190 

4.712 

2.362 

2.403 

-- 

2.347 

3.306 

. 986 

2.438 

200 

4.608 

2.288 

2.174 


2,315 

3.267 

.988 

2.614 

210 

4.551 

2.271 

2.120 

-- 

2.278 

3.213 

. 989 

2.790 

220 

4.521 

2.298 

2.207 

! 

2.221 

3.136 

. 987 

2 . 965 

230 

4.512 

2.359 

2.393 


2 . 151 

3.040 

. 981 

3.139 

240 

5.734 

3.575 

4.934 

— 

2.156 

3.048 

0.983 

1.810 

170 

5.217 

3.100 

4.089 

-- 

2.114 

2.989 

. 983 

1.988 

180 

4.850 

2.750 

3.374 

-- 

2.098 

2.966 

. 984 

2.106 

190 

4.599 

2.506 

2.797 

-- 

2.091 

2. 956 

. 985 

2.343 

200 

4.435 

2.352 

2.373 

-- 

2.080 

2.941 

. 986 

2.521 

210 

4.333 

2.270 

2.117 

-- 

2.060 

2.914 

. 988 

2.697 

220 

4.27 9 

2.245 

2.033 

-- 

2.031 

2.872 

. 990 

2.673 

230 

4.254 

2.265 

2.099 

-- 

1.987 

2.809 

. 988 

3.048 

240 

4 . 254 

2.319 

2.273 


1.932 

2.730 

.983 

3.222 

250 

5.490 

3.544 

4.881 

-- 

1.944 

2.748 

0.983 

1.894 

180 

4.985 

3.086 

4.061 

-- 

1.897 

2.678 

. 983 

2.072 

190 

4.623 

2.744 

3.361 

-- 

1.877 

2.649 

.984 

2.250 

200 

4.374 

2.502 

2.787 

-- 

1.869 

2.636 

.985 

2.427 

210 

4.209 

2.345 

2.353 

-- 

1.861 

2.625 

. 986 

2.605 

220 

4.108 

2.258 

2.078 

-- 

1.847 

2 . 604 

.988 

2.781 

230 

4.055 

2.227 

1.970 

__ 

1.825 

2.572 

. 990 

2 . 957 

240 

4.034 

2.240 

2.016 

-- 

1.792 

2.521 

.989 

3.132 

250 

4.041 

2.289 

2.179 


1.749 

2.456 

.984 

3.306 

260 

5.298 

3.508 

4.820 

__ 

1.788 

2.515 

0 . 983 

1.979 

190 

4.801 

3.068 ! 

4.030 

-- 

1.731 

2.427 

. 983 

2.157 

200 

4.443 

2.736 

3 344 

; 

1 . 704 

2.385 

. 984 

2.335 

210 

4.193 

2.499 j 

2.777 

-- 

1.692 

2.366 

. 985 

2.513 

220 

4.028 

2.341 ! 

2.339 

i 

1.684 

2.354 

. 986 

2.690 

230 

3.926 

2 . 250 

2.050 


1.673 

2.336 

. 988 

2.866 

240 

3.874 

2.214 

1.923 


1.65 7 

. 2.311 

. 990 

3.042 

250 

3.857 

2.222 

1. 953 


1.632 

2.271 

. 98 9 

3.217 

260 

3.872 

2.267 

2.108 

-- 

1.600 

2.218 

.985 

3.391 

270 

5.799 

4.003 

5.621 

__ 

1.796 

2.527 

0.984 

1.888 

190 

5. 148 

3.469 

4.754 

-- 

1.677 

2.343 

.983 

2.066 

200 

4.659 

3.048 

3.988 

__ 

1.610 

2.233 

.983 

2.244 

210 

4 . 303 

2.726 

3.321 

-- 

1.575 

2.175 

. 984 

2.422 

220 

4.053 

2.493 

2.763 

-- 

1.557 

2.146 

.985 

2.599 

230 

3.885 

2.336 

2 . 325 

-- 

1.547 

2.128 

.986 

2.776 

240 

3.782 

2.244 

2.029 

-- 

1.536 

2.111 

. 988 

2.953 

250 

3.730 

2.204 

1.889 

l 

1 . 524 

2.090 

. 990 

3.128 

260 

3.717 

2.209 

1 . 907 


1.505 

1 2.058 

i 

.990 

3.303 

270 

5.672 

3 . 935 

5.514 

-- 

1.735 

2 434 

0.984 

1.976 ! 

200 

5.034 

3.426 

4.679 

-- 

; 1 . 606 

2 . ?27 

.983 

2.154 

210 

4.552 

3.023 

3.939 

-- 

i 1.527 

2.095 

.983 

2.332 

220 

4.197 

2.713 

3 . 29 ? 

-- 

1.482 

2.010 

.904 

2 . 509 

230 

3.946 

2.486 

2.744 

-- 

1.458 

1 . 975 

. 985 

2.687 

240 

3.777 

2.331 

2.510 

__ 

1.444 

1 . 94 9 

. 986 

2.864 

250 

3.673 

2 238 

2.010 

-- 

1 . 4 33 

1 . 930 

.988 

3 . 04 1 

260 

3.621 

2.198 

1.864 


1.4 23 

1.911 

. 990 

3.216 

270 

3.612 

2.201 

1.075 

-- 

1.411 

1.087 

.990 

3.391 

280 




TABLE II. 


CHARACTERISTICS OP OPTIMAL MARS -TO -EARTH TRAJECTORIES FOR 1980 


Julian date 
at Mars 
departure 

Transit 
time, T rf _ e , 
days 

Total velocity 
increment, 

£! AV i- 
miles/sec 

Mars de 

parture 

Kick, 

Earth arrival 

Minimum 

Central 
angle, (3, 
rad 

av 4, 

miles/sec 

V co 4 , 

miles/ sec 

av 5 , 

miles/sec 

AV 6, 

miles/sec 

V “6, 

miles/sec 

radius, 

^mim 

astronomical 

units 

2444420 

220 

9.771 

3.611 

4.862 

1.566 

4.595 

6.498 

0.548 

4.912 


230 

9.202 

3.414 

4.631 

1.767 

4.022 

5.650 

.546 

5.089 


240 

8.745 

3.189 

4.364 

1.802 

3.753 

5.225 

.540 

5 . PfiS 


250 

8.369 

3.074 

4.225 

1.924 

3.371 

4.583 

.539 

5.440 


260 

8.081 

2.968 

4.095 

2.001 

3.113 

4.113 

.539 

5.615 


270 

7.860 

2.915 

4.030 

2.118 

2.826 

3.538 

.543 

5.788 


280 

7.716 

2.830 

3.925 

2.175 

2.712 

3.290 

.542 

5 . 961 


290 

7.620 

2. 724 

3.793 

2.298 

2.598 

3.025 

.552 

6.132 


300 

7.597 

2.696 

3.758 

2.358 

2.543 

2.890 

.550 

6.303 

2444430 

210 

10.424 

4.004 

5 . 315 

i.526 

4.894 

6. 918 

0.542 

4.828 


220 

9.749 

3.650 

4.907 

1.617 

4.482 

6.336 

.531 

5.005 


230 

9.190 

3.490 

4.721 

1.829 

3.870 

5.413 

.531 

5.181 


240 

8.744 

3.263 

4.453 

1,839 

3.642 

5.044 

.527 

5.357 


250 

8.387 

3.144 

4.310 

1.936 

3.307 

4.469 

.526 

5.531 


260 

8.112 

3.077 

4.228 

2.056 

2.979 

3.853 

.529 

5.704 


270 

7.922 

2. 982 

4.113 

2.115 

2.825 

3.537 

.529 

5.877 


280 

7.880 

2.731 

3.802 

2.383 

2.766 

3.410 

.553 

6.048 


290 

7.738 

2.831 

3.926 

2.304 

2.603 

3.038 

.536 

6.219 

2444440 

210 

10.387 

3.972 

5.278 

1.509 

4.905 

6.933 

0.521 

4.920 


220 

9.735 

3.650 

4.908 

1.572 

4.514 

6.381 

.515 

5 . 096 


230 

9.182 

3.570 

4.814 

1.870 

3.742 

5.208 

.518 

5.271 


240 

8.746 

3.395 

4.609 

1.938 

3.413 

4.656 

.516 

5 . 446 


250 

8.416 

3.220 

4.401 

1.944 

3.252 

4.370 

.512 

5.619 


260 

8.164 

3.149 

4.316 

2.053 

2.961 

3.817 

.515 

5.792 


270 

7.999 

3.054 

4.201 

2.122 

2.823 

3.532 

.516 

5.963 


280 

7.901 

2.992 

4.126 

2.195 

2.713 

3.293 

.516 

6 . 134 


290 

7.892 

2.916 

4.032 

2.352 

2.625 

3.089 

.526 

6.303 

2444450 

200 

11.136 

4.301 

5.650 

1.309 

5.526 

7.770 

0.508 

4.833 


210 

10.353 

4.020 

5.333 

1.561 

4.772 

6.748 

.505 

5.010 


220 

9.708 

3.732 

5.003 

1.643 

4.333 

6.118 

.500 

5.185 


230 

9.182 

3.555 

4.797 

1.768 

3.858 

5.394 

.499 

5.35 9 


240 

8.766 

3.425 

4.645 

1.877 

3.464 

4.744 

.499 

5.533 


250 

8.449 

3.337 

4.541 

1.991 

3.120 

4.127 

.502 

5. 705 


260 

; 8.231 

3.226 

4.409 

2.048 

2.956 

3.808 

.502 

5.877 


! 270 

8 . 108 

3.177 

4.349 

2.095 

2.836 

3.560 

.500 

6.047 


280 

8 . 033 

3.058 

4.205 

2.293 

2.683 

3.223 

.515 

6.217 

2444460 

200 

! 11.093 

4.391 

5.752 

1.434 

5.269 

7.428 

0.494 

4.922 


210 

10.326 

4.030 

5.344 

1.527 

4.769 

6.744 

.488 

5.097 


220 

9.689 

3.821 

5.105 

1. 692 

4.177 

5.885 

.486 

5.271 


230 

9.186 

3.638 

4.893 

1.782 

3.766 

5.247 

.486 

5.445 


240 

8.792 

3.514 

4.749 

1.885 

3.393 

4.621 

.487 

5.617 


250 

8.505 

3.416 

4.633 

1.979 

3.110 

4.108 

.489 

5.789 


260 

8.316 

3.309 

4.507 

2.042 

2.965 

3.825 

.487 

5. 959 


270 

8.204 

3.245 

4.431 

2.148 

2.811 

3.507 

.491 

6.129 


280 

8.197 

3.103 

4.260 

2.405 

2.608 

3.236 

.510 

6.298 

2444470 

200 

11.057 

4.505 

5.879 

1.543 

5.009 

7.077 

0.482 

5.008 


210 

10.298 

4.108 

5.433 

1.509 

4.621 

6.535 

.474 

5.182 


220 

9.680 

3.964 

5.269 

1.783 

3.933 

5.511 

.47 7 

5.356 


230 

9.199 

3.733 

5.003 

1.795 

3.672 

5.093 

.473 

5.528 


240 

l 8,833 

3.608 

4.859 

1.889 

3.336 

4.522 

.475 

5.699 


250 

8.576 

3.500 

4.732 

1.969 

3.107 

4.102 

.475 

5.870 


260 

8.429 

3.436 

4.658 

2.026 

2.967 

3.829 

.470 

6.040 


270 

8.360 

3.311 

4.510 

2.279 

2.772 

3.423 

.493 

6.208 


280 

8.388 

3.121 

4.282 

2.565 

2. 702 

3.267 

.506 

6.377 

2444480 

200 

11.016 

4.495 

5.868 

1.494 

5.027 

7.100 

0.464 

5.091 


210 

10.275 

4.213 

5.551 

1.624 

4.439 

6.272 

.462 

5.265 


220 

9.684 

4.055 

5.373 

1.788 

3.841 

5.366 

.464 

5.437 


230 

9.227 

3.874 

5.166 

1.848 

3.505 

4.815 : 

.464 

5.609 


240 

8.890 

3.712 

4.980 

1.893 

3.285 

4,430 

.463 

5.779 


250 

8.669 

3.592 

4.840 

1.957 

3.120 

4.127 

.461 

5.949 


260 

8.538 

3.493 

4.724 

2.157 

2.889 

3.669 

.476 

6.118 


270 

8.546 

3.331 

4.534 

2.440 

2.774 

3.428 

.491 

6.286 


280 

8.588 

3.080 

4.233 

2.787 

2.721 

3.309 

.501 

6.453 


2444490 


200 

210 

220 

230 

240 

250 

260 

270 

280 


10.988 

4.563 

10.271 

4.268 

9.703 

4.091 

9.271 

3. 963 

8.967 

3.805 

8.786 

3.706 

8.761 

3.567 

8.773 

3.326 

8.807 

2. 937 


5.943 

1.508 

5.614 

1.590 

5.413 

1.713 

5.268 

1.831 

5.071 

1.877 

4.972 

1.942 

4.811 

2.349 

4.527 

2.665 

4. 057 

3.092 


4.918 

6.951 

4.413 

6.235 

3.897 

5.454 

3.477 

4.767 

3.284 

4.428 

3.138 

4.161 

2.845 

3.580 

2.782 

3.445 

2.779 

3.437 


0.450 

5.172 

.448 

5.345 

.447 

5,516 

.451 

5.687 

.450 

5.856 

.446 

6.025 

.480 

6.193 

.489 

6.361 

.493 

6.528 


4.981 

4.641 

4.337 

4.189 

4.056 

3.905 

3.804 

3.591 

3.263 


6.405 

6.031 

5.691 

5.524 

5.362 

5.201 

5.085 

4.839 

4.452 


1.374 

1.513 

1.558 

1.703 

1.796 

1.864 

2.119 

2.501 

2.873 


5.483 

4.817 

4.387 

3.841 

3.491 

3.296 

3.003 

2.867 

2.823 


7.713 

6.810 

6.197 

5.366 

4.791 

4.450 

3.901 

3.624 

3.533 


0.439 

.437 

.436 

.436 

.438 

.436 

.455 

.474 

.481 


5.078 

5.251 

5.422 

5.593 

5.762 

5.931 

6.099 

6.268 

6.434 



Total impulsive velocity increment, Total impulsive velocity increment, 
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Figure 8. - Comparison of coplanar and three-dimensional results 
for Earth -Mars trips. Earth -Mars travel time, 200 days. 


normal to the ecliptic in order to 
include the Earth, Sun, and Mars ; 
the associated Y, AV is unreal- 
istically high. The optimal 
three-dimensional trajectory for 
this case is a broken-plane trans- 
fer as described in reference 5 
and indicated by the dotted curve 
on figure 8. It may be seen from 
this comparison that the present 
coplanar results are actually a 
better approximation to optimal 
three-dimensional transfers than 
are the single -plane transfers of 
reference 2. 

Characteristics of double - 
conic transfers . - The inward 
transfers shown in table II in- 
volve long travel angles; they are 
double-conic trajectories in which 
the midcourse impulse was chosen 



fers result in a sizable reduction 

in ^ 4 AV as compared to s ingle - 

conic trajectories. This is il- 
lustrated in figure 9, where 
Y, AV is plotted against Mars 
departure date (and travel angle) 
for 260-day single-conic (dashed 
line) and double-conic inward 
transfers (solid curve) . As was 
the case with the circular - 
coplanar data used in illustrating 
the search procedure, the double- 
conic transfers become advanta- 
geous for travel angles greater 
than about 4 radians; the corre- 
sponding departure date is 2444220. 
The advantage increases rapidly 
for larger travel angles (later 
departure dates). For instance, 

a double-conic 260-day trip departing on 2444460 ((3 = 5.96 rad) requires a 
Y AV of only 8.3 miles per second, while the equivalent single-conic requires 
nearly 20 miles per second. 


to minimize 




These trans- 


2444 160 260 360 

Julian date at Mars departure 


460 


Figure 9. - Characteristics of double-conic Mars-Earth transfers 
in 1980. Mars-Earth travel time, 260 days. 
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Julian date at Mars arrival 


Figure 10. - Mars round-trip optimization. Total trip time, 500 days; 
stay time, 40 days; low circular parking orbits. 

Round Trips 

The optimal one-way transfer data presented in tables I and II may be used 
for the construction of round trips. All that remains is to select the best 
distribution of travel time and angle between the outward and return legs, for 
each desired trip time and stay time at Mars. The selection of the travel time 
and angle distribution is equivalent to the selection of the date of Mars en- 
counter and the outward travel time. Consider, for instance, a 500-day round 
trip with a 40-day stay at Mars. Choosing an arrival date = 2444420 and 

AT 0-C r = 200 days, it may be seen in table I that S q AV = 4.599 miles per sec- 
ond for the outward transfer (assuming propulsive braking at Mars). The return 
transfer thus departs Mars on 2444460 and has a travel time of 260 days) from 

table II, the corresponding ^ ^ AV = 8.316 miles per second, with propulsive 
braking at Earth. Hence, for this particular round trip, the Sq ^ ~ 12.915 
miles per second. By repeating this procedure, the most favorable arrival date 
and travel time may be found as illustrated in figure 10. The ^ AV for all 
propulsive single-conic (dashed curves) and double-conic (solid curves) round- 
trip durations of 500 days with 40-day stay is plotted against Mars arrival date. 
In each case, the minimum value is denoted by the circle symbols. The full AV 
savings of 3.2 miles per second is developed at an arrival data and travel times 
considerably different from those which were best for single-conic transfers (a 
modest saving, the square symbol, is still available when double -conic transfers 
are substituted directly into a minimum AV single -conic round trip without 

reoptimizing the dates and travel times) . The overall effect of the reoptimi- 
zation is to allocate a larger travel time and angle to the long-angle return 
transfer, thus taking full advantage of the favorable low ^ AV at long 
travel angle characteristic of double-conic transfers. 

The minimum ^ AV single- and double-conic trajectories in figure 10 
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TABLE III. - COMPARISON OF OPTIMAL SINGLE- AND DOUBLE-CONIC ROUND TRIPS IN 1980 
[Total trip time, 500 days; stay time, 40 days.] 


Type of 
trajectory 

Minimum 

Outward 

transfer 


Return transfer 

bdJL 

velocity 

increment, 

Ei AV, 

Earth- 

Mars 

travel 

time, 

AT ®^, 

days 

Earth- 

Mars 

travel 

angle, 

Pe-cr, 

rad 

Impulsive velocity 
increment, 
miles/sec 

Mars- 

Earth 

travel 

time, 

AT C^B, 

days 

Mars- 

Earth 

travel 

angle, 

rad 

Impulsive velocity 
increment , 
miles/sec 


miles/sec 

At 

Earth, 

AV i 

At 

kick 

point, 

av 2 

At 

Mars, 

AV 3 

At 

Mars, 

av 4 

At 

kick 

point, 

av 5 

At 

Earth, 

AV 6 

Single-conic 

16.3 

260 

3.13 

2.20 

0 

1.50 

200 

5.17 

4.25 

0 

8.35 

Double -conic 

12.9 

200 

2.34 

2.40 

0 

2.30 

260 

5.96 

3.25 

2.0 

2.95 
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Figure 11. - Comparison of single- and double-conic round trips to Mars. All 
propulsive braking (1980). 


are further compared in table III. It is of interest to note that, with single- 
conic transfers, the individual velocity increments vary from 1.5 to 8.35 miles 
per second, while with double-conics the variation is only between 2 and 3.25 
miles per second. It may be inferred from the AVg data (Earth return) that 
an extremely high (about 70 000 ft/sec) Earth aerobraking capability would be 
required to make the single-conic trajectory competitive with the double-conic. 

Effect of total trip time . - In the preceding example, the use of double- 
conic transfers was seen to produce a large reduction in 2^ for a 500-day 

trip time. It is also of interest to consider how the available savings varies 
as a function of the total trip time. In figure 11, the ^ AV is plotted 
against total trip time for all -propulsive, single- and double-conic trajecto- 
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ries in 1980. A stay time of 40 
days was chosen to be consistent 
with existing single-conic round- 
trip data (ref. 2). It will 
be noted that the single-conic 
trajectories (dashed lines) are 
represented by three distinct 
curves. The first curve (at 
shortest trip times) corresponds 
to fast trips in which the space- 
craft describes a total central 
angle equal to the motion of the 
Earth (N - 0; cf. eq. (ll)). The 
second curve, also with 40-day 
stays, involves longer trip times 
and the spacecraft travels exactly 
one revolution less than the 
Earth (N « -l) . The transition 
region between these two curves 
involves central angles that are either very long or very short. The last curve 
represents near-Hohmann trips with N = -1 and stay times of 450 days. It is 
evident from the figure that double-conic trajectories (represented by the solid 
line) yield significantly lower XI AV than the single-conic ones in the 
"fast" (N = 0) trip region. The saving becomes appreciable at about 380 days 
and increases to 6.5 miles per second out of 19 miles per second at 640 days. 

A new but very shallow minimum may be seen at 560 days, where XI AV - 12.3 
miles per second, a saving of 3.8 miles per second over the single-conic tra- 
jectory minimum. At trip times approaching 700 days, single-conic trajectories 
represented by the second dashed curve (N = -l) become competitive with the 
double-conic trajectories. The double conics are not advantageous for trip 
times greater than about 700 days. The lowest X) AV of all corresponds to 
Hohmann-type round trips with long trip times (850 to 1000 days) and long stay 
times at Mars (e.g., 450 days). 

Effect of stay time at Mars . - In figure 12, the ^ AV is plotted against 
the stay time at Mars for 420-, 500-, 600-, and 700-day round trips. The solid 
curves represent double-conic trajectories while single-conics are denoted by 
the dashed curves. Here it is evident that the double-conic trajectories re- 
duce JI AV at a given trip time over a wide range of stay times. The reduc- 
tion, however, is somewhat less pronounced at the longer stay times. It is of 
interest to note that if double-conic trajectories are used, a 500-day trip with 
a 130-day stay has the same XI AV as a 500-day single-conic trip with a stay 
time of only 40 days . 

The dot-dashed envelope curve on the lower part of the figure represents 
the lowest AV obtainable with given stay times using double-conic trajec- 

tories in the "fast” (N = 0) trip region. An additional parameter, the trip 
time, was left open in obtaining this curve. The optimal values of T^ ^ cor- 
responding to these trip times, as indicated by the markers along the envelope 
curve, range from 540 days (T_ = 0) to 640 days (T = 140) . There is clearly 

very little incentive for considering trip times longer than about 550 days 



Figure 12. - Effect of stay time at Mars on minimum total velocity 
increment. Optimum single- and double-conic trips in 1980. 
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Delay in Mars departure, days 

Figure 13. - Effect of delayed departure from Mars on minimum total velocity 
increment. Optimum single- and double-conic trips in 1980. Total trip 
time, 500 days; stay time, 40 days; Earth atmospheric entry velocity, 

V 0e < 45 000 feet per second. 


unless very long stay times are required. 

Delayed departure from Mars . - The trajectories discussed up to this point 
are characterized by the use of advance planning to select the most favorable 
transit time and central angle distribution for each chosen value of ^ and 

T s . A realistic flight plan, however, should be flexible enough to permit a 
safe return in spite of a minor mishap. Delayed departure from Mars will be 
discussed here as a typical example of the many types of safe-return problems 
that might be considered. This could become necessary for many plausible rea- 
sons, which include (l) rendezvous difficulty in returning a landing partly to 
the mother ship and (2) a minor propulsion system malfunction during nuclear 
engine startup. Either of these problems could require a considerable time de- 
lay for correction or repair. 

The problem of delayed departure is similar to an extension of stay time 
except that only the inbound transfer can be reoptimized in response to the 
changed trajectory requirements. Figure 13 compares the AV increase for 

single- and double-conic trajectories for delays in departure ranging up to 
70 days. The comparison is based on planned minimum AV single- and double- 
conic trips in 1980 with a 500-day duration and a 40-day stay time at Mars. Two 
cases were considered for each type of trajectory: First it was assumed that 

the round trip, including the delay, must be completed within the original 500- 
day period, and second that the round-trip time can be extended by as much as 
60 days (to 560 days). With no time extension, both single- and double-conic 
trajectories require increased AV for late departure, but the increase is 
considerably reduced if double-conic trajectories are used. If the trip time 
can be extended to 560 days, the increase is smaller for both types, but again 
the increase for the double-conic trajectory is considerably less than it is 
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Figure 14. - Effect of atmospheric braking at Earth. Earth entry velocity, 
V 0e < 52 000 feet per second on round trips to Mars. 


for single conics. This comparison shows that the ^ AV increase for depar- 
ture delays can he cut by a factor of about 2 by using double-conic trajecto- 
ries. This comparison is for trajectories which use propulsive braking at Earth, 
or atmospheric braking from velocities up to 45 OOO feet per second. The dif- 
ference is smaller for higher entry velocities. 

Atmospheric braking . - It was shown above that double-conic trajectories 
can be used to produce large reductions in the all -propulsive ^ AV at trip 
and stay times of interest and to mitigate the effect of unexpected delay in 
Mars departure. It is well known that atmospheric braking can also be used to 
reduce £ AV. This section will consider the propulsive velocity increments 
and atmospheric entry velocities that occur when double-conic trajectories are 
used in combination with atmospheric braking at Earth return. 

The combination of atmospheric braking with double-conic trajectories may 
be assessed according to two criteria, that is, the AV savings compared to 

the all -propulsive double-conic case and compared to the single-conic trajectory 
with equal atmospheric braking. The results of these comparisons depend entirely 
on the allowable entry velocities . Assuming 37 OOO feet per second as a limit 
(Apollo technology), single- and double-conic trajectories would benefit equally; 
the curves that were shown in figure 11 would simply be shifted downward by 
about 2 miles per second without affecting the relative advantage of double - 
conic trajectories. If the entry velocity is unrestricted, however, the double- 
conic trajectories are not appreciably better than single conics. In this case 
the search and optimization procedure described previously yields results ap- 
proaching single-conic transfers. The associated entry velocities, however, are 
extremely high (e.g., about 70 000 fps for the 500-day trip). Assuming an in- 
termediate limit, V <52 000 feet per second, leads to the results shown in 
figure 14. Here AV is plotted against T tot as it was in figure 10; com- 
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Figure 15. - Effect of allowable Earth atmospheric entry velocity on minimum 
total velocity increment for 500-day round trips to Mars in 1980. Stay 
time, 40 days. 


parison of these two figures indicates that the AV reduction due to double- 

conic trajectories is available over the same wide range of trip times as in the 
all -propulsive case. The minimum E AV for the double-conic is 2.5 miles per 
second below the minimum available from a single -conic trajectory with the same 
entry velocity limit. This is also a saving of 3.5 miles per second as compared 
to the all -propulsive double-conic case. 

The variation of ^ AV with allowable Earth atmospheric entry velocity 
is further illustrated in figure 15, where the 500-day trip with a 40-day stay 
is used as an example. The E AV is plotted against V^ using dashed and 

solid curves to represent single- and double-conic trajectories, respectively. 

The maximum advantage (3.2 miles/sec) that occurs in the all -propulsive braking 

case is maintained up to an entry velocity of about 45 000 feet per second, and 

thereafter decreases as V is further increased. A retromaneuver is needed 

e® 

to reach entry velocities below 45 000 feet per second. The orbit elements of 
the Mars -Earth transfer are approximately constant for entry velocities from 
25 000 to 45 000 feet per second; the reduction in E ^ as ^e® increases 

in this range is simply due to the proportional decrease in the retromaneuver 
AV. If V e0 > 45 000 feet per second, no retromaneuver is used; instead, the 
transfer orbit elements are reoptimized to take full advantage of the higher 
V e® limit * 

From this example it may be observed that double-conic trajectories can be 
used in combination with atmospheric braking to produce attractive reductions in 
E AV without requiring excessive entry velocities. The advantage of double- 
conic over single-conic trajectories tends to decrease as the allowable Earth 
atmospheric entry velocity is increased. 
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Figure 16. - Comparison of single- and double-conic round trips to Mars. All 
propulsive braking (1971k 


Effect of synodic period . - Because of the eccentricity of Mars's orbit, 
its radius at opposition varies from 1.38 astronomical units in 1971 to 1.66 as- 
tronomical units in 1980. Optimal trips in the latter period involve greater 
radial and angular travel and thus tend to require a higher 23 than trips 

in other oppositions. Conversely, trips in 1971 are the least difficult. Trips 
in the difficult 1980 period were shown in the preceding figures. Figure 16 il- 
lustrates the variation of 23 AV with T ^ ^ for all -propulsive trips in 1971 

using solid lines for double-conic trajectories and dashed lines for single con- 
ics. Again, the double-conic trajectories yield appreciable reductions of 
Y. AV for trip times of 360 to 720 days. The reductions, however, are less 
pronounced than those found in 1980. The minimum 23 AV is decreased by 
1.6 miles per second in this period (compared to 3.8 miles/sec in 1980). 

As was the case in 1980, there is little £ AV advantage for double-conic tra- 
jectories if unrestricted Earth aerobraking capability is available. 

It is pointed out in appendix B that the advantage of a double-conic over a 
single-conic trajectory increases as the angular and radial travel increase. 

Thus, the same factors that make 1971 an "easy" year also account for the less 
advantageous showing of double-conic trajectories. 

Five -hundred-day single- and double-conic round trips are compared in fig- 
ure 17 over a range of synodic periods . The AV sum for the minimum 23 AV 
all-propulsive trips within each synodic period is plotted against opposition 
year from 1971 to 1986. The double-conic trajectories produce a ^ AV reduc- 
tion in all synodic periods, ranging from 1.2 miles per second in 1971 and 1986, 
to 3.2 miles per second in 1980. The net result is a pronounced flattening of 
the 23 AV requirement as a function of synodic period in addition to the gen- 
eral reduction in level. 

Effect of trajectory profile . - It was illustrated in table III that the 
optimum travel angle and travel time distribution for fast (N « 0) Mars round 
trips are not necessarily symmetrical, and also that the assymetry is increased 
when the trajectory is reoptimized for double-conic transfers. The trips that 
were discussed for 1980 consisted of short-angle outward transfers and long- 
angle returns. These are termed "short-long" trips. Another class of round 
trips, "long-short," also exists, in which the angular relation is reversed. 
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Figure 17. - Comparison of optimum single- and double-conic round trips to 
Mars in different synodic periods. Total trip time, 500 days; stay time, 

40 days. 


It is evident from figure 17 that the choice of trajectory profile has a 
significant effect on 53 AV in most synodic periods. For both single- and 
double-conic trajectories, the long-short profile yields minimum all -propulsive 
53 AV from 1971 to 1980, while the short-long profile is advantageous between 
1980 and 1986. (As shown in appendix B, the more difficult profile in a given 
period is the one which requires the greater radial and angular travel.) 

For the trip time chosen in this example, the double-conic trajectories 
possess distinct local -minimum long-short and short-long profiles in every syn- 
odic period, (in 1980, for instance, the local minimum for the short-long pro- 
file is illustrated by the lower circle symbol on fig. 10.) In most years only 
one trajectory profile yields a local minimum 53 AV for single-conic trips. 

These facts have several implications. In a given synodic period, double- 
conic trajectories not only reduce the ^3 AV , but also reduce the variations 
between the individual AV's. The ^ AV penalty in designing for the worst 
synodic period compared to the best, is only 2 miles per second, if double- 
conics are used. Moreover, the individual AV distribution, being relatively 
uniform, would not be significantly affected when the transition from the long- 
short to the short-long profile is made. 

By contrast, with single-conics, the 53 penalty for designing for the 

worst period is over 4 miles per second. In addition, the individual AV dis- 
tribution, being highly nonuniform (cf. table III), would be drastically modi- 
fied in making the transition from the long-short to the short-long trajectory 
profile. 

These considerations suggest that by using double-conic trajectories, a 
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standardized spacecraft could be designed which could accomplish a mission in 
any synodic period, or a sequence of missions covering many synodic periods. 

Comparison with Venus swingby trajectories . - Another efficient new class 
of Mars round-trip trajectories is the Venus swingby described in references 9 
and 10. These trajectories incorporate a close passage by Venus during the 
long-angle leg of a round trip. They resemble the present trajectories in that 
low ^2 AV long -angle transfers consist of two conic arcs; the gravitational 
force acting during the Venus passage replaces the midcourse propulsion as the 

TABLE IV. - COMPARISON OF DOUBLE-CONIC AND VENUS SWINGBY TRAJECTORIES 

FOR ROUND TRIPS TO MARS 

[Propulsive braking ; low circular parking orbits; 
stay time, 10 days; optimum total trip time.] 


Launch 

year 

Type of 
trajectory 

Minimum 

total 

velocity 

increment, 

miles/sec 

Optimum 

total 

trip 

time, 

T tot' 

days 

Velocity increments 

Outbound transfer 

Return transfer 

AV i 

(date, 

244-) 

AV 2 

(date, 

244-) 

AV 3 

(date, 

244-) 

AV 4 

(date, 

244-) 

AV 5 

(date, 

244-) 

AV 6 

(date, 

244-) 

a 1975 

Single- 

conic 

16.37 

430 

2.58 

(2670) 


2.68 

(2840) 

3.14 

(2850) 


7.97 

(3100) 

b 1975 

Double - 
conic 

10.81 

510 

2.50 

(2366) 

1.60 

2.71 

(2642) 

1.87 

(2652) 


2.13 

(2867) 

a 1975 

Venus 

swingby 

11.18 

518 

2.61 

(2670) 


3.35 

(2820) 

2.24 

(2830) 

Pass 

Venus 

(3050) 

2.98 

(3188) 

a 1978 

Single- 

conic 

16.80 

430 

2.51 

(3440) 


2.91 

(3620) 

3.15 

(3630) 


8.23 

(3870) 

b 1978 

Double - 
conic 

11.49 

517 

2.53 

(3109) 

2.06 

2.77 

(3397) 

1.86 

(3407) 


2.27 

(3625) 

a 1978 

Venus 

swingby 

13.05 

699 

3.23 

(3081) 

Pass 

Venus 

(3319) 

3.23 

(3458) 

1.87 

(3468) 


4.72 

(3780) 

a 1980 

S ingle - 
conic 

15.77 

430 

2.42 

(4190) 


3.18 

(4380) 

2.96 

(4390) 


7.21 

(4620) 

b 1980 

Double - 
conic 

11.88 

550 

2.22 

(4180) 


1.63 

(4440) 

3.06 

(4450) 

2.29 

2.68 

(4730) 

a 1980 

Venus 

swingby 

10.40 

540 

2.85 

(3840) 

Pass 

Venus 

(3996) 

3.45 

(4120) 

1.45 

(4120) 


2.65 

(4380) 


^ef. 10; minimum weight. 

^Present analysis; minimum T] ^ AV. 
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means of switching from one arc to the next. 

The Venus swingby trajectories presented in reference 10 are not directly 
comparable to the present results because they were obtained on the bais of min- 
imum gross weight for a particular space vehicle configuration and for a mission 
profile (atmospheric braking at Mars) that was not studied herein. A limited 
comparison (table IV) nevertheless yields some interesting similarities and con- 
trasts: 

(1) In comparison with single-conic trajectories, both double-conic and 
Venus swingby trajectories offer significant 53 AV savings together with mod- 
erate Earth-approach velocities. Only one of the double-conic or swingby tra- 
jectories shown in table IV (the 1978 swingby) would benefit from an Earth aero- 
braking capability greater than about 42 000 feet per second. The single conics 
would need 65 000 to 75 000 feet per second of aerobraking capability to be com- 
petitive. 

(2) In 1975, the optimum double-conic and Venus swingby trajectories will 
require comparable values of ^ ® AV and trip time. 

(3) In 1978, the double-conic trajectory appears to have a clear advantage 
from the standpoints of both £ p AV and trip time. 

(4) In 1980, by contrast, the Venus swingby offers a decidedly lower value 

Z 6 

^ AV for comparable trip times. 

(5) It may be noted that there is a considerable difference in the Mars- 
arrival dates (and hence, the Martian season) for double-conic and swingby tra- 
jectories. This may be significant from the standpoint of mission objectives. 

Although 1980 is the most difficult period for double-conic trajectories, 
it is especially favorable for swingbys. In such years, Venus is in the most 
appropriate position relative to Earth and Mars. It is then possible to obtain 
many of the characteristics of double-conic trajectories without midcourse pro- 
pulsion. In other years, however, the elimination of midcourse propulsion is 
offset by less favorable configurations of Venus. 

On the other hand, double-conic trajectories owe their efficiency to a 
careful optimization of the location, magnitude, and direction of the midcourse 
AV. This process does not depend on the configuration of a third planet and 
therefore yields consistent results in all synodic periods. 

These examples indicate that double-conic and Venus swingby trajectories 
are comparable techniques whose relative merits depend strongly on the synodic 
period and possibly on mission objectives. Further comparisons between them 
will depend on accounting for both mission factors and trajectory character- 
istics in a consistent manner. Finally, the possibility cannot be excluded that 
hybrid trajectories, which would incorporate both Venus swingby and double-conic 
portions, may yield a lower 53 AV than either double conics or Venus swingbys 
alone . 
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CONCLUSIONS 


An efficient class of high-thrust interplanetary trajectories has been an- 
alyzed in which a midcourse impulse is used to reduce the total velocity re- 
quirement. Each one-way plane-to-planet transfer consists of two heliocentric 
conic sections and requires three impulses rather than the customary two. 

For one-way trips to Mars, these double-conic trajectories yield lower to- 
tal velocity increments than conventional single-conic trajectories if the 
heliocentric travel angle is greater than about 4 radians. The improvement is 
most noticeable at long transit times and is large enough to produce new rel- 
ative minima in the long -angle region for transit times greater than 350 days. 
None of the new time-constrained minima, however, are lower than that resulting 
from the classical single-conic Hohmann transfer. 

Although Hohmann-type round trips of about 3-year duration represent an 
absolute minimum AY, the double-conic trajectories can be used to produce 

notable reductions in ^ AV for round trips of 1 to 2 years. The reductions 
are largest for all propulsive trips but are still significant when atmospheric 
braking is used at Earth return. The greatest savings occur in the most diffi- 
cult synodic periods; consequently, the variation of ^ AV with synodic peri- 
od is decreased by a factor of about 2. The variation of the individual veloc- 
ity increments between long-short and short-long trips is also reduced by a 
large amount. These properties of double-conic trajectories may imply gross- 
weight reductions for an interesting class of trips and enhance the possibility 
of using a standard spacecraft for missions in many synodic periods. 

A comprehensive mission study is required to more precisely evaluate the 
advantages and applicability of double-conic trajectories as compared with Venus 
swingby and conventional single-conic trajectories for specific Mars missions. 
Additional effort is indicated to investigate the utility of double-conic tra- 
jectories for missions to Venus and the major planets, and to determine whether 
two or more midcourse impulses could be used to obtain still greater reductions 
in £ AV. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 28, 1965. 
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APPENDIX A 


COMPUTATION OF PLANET POSITION AND VELOCITY 


Elliptic Planet Orbits 


The elements of a transfer trajectory, and the associated velocity incre- 
ments, may be computed in terms of the travel time if the position and velocity 
vectors at each planet terminal are known. These quantities may in turn be de- 
rived from the date of encounter at each planet terminal by using the planet- 
orbit data shown in table V (obtained from ref. 3) . 

Only one more parameter, the true anomaly 0 (d), is required to compute the 
planet position and velocity at the encounter date. The elapsed time since the 
last perihelion passage AT is 

AT a D - D + NT (Al) 

p pr ' ' 


where N( = 0, ±1, ±2, . . . ) is chosen to obtain the smallest positive value of 
ATp. An assumed or estimated value of 0 (d), inserted in the time equation (5), 
yields a trial value that may be compared to equation (Al) . The error is then 
used to obtain a better estimate of 0 (D); this iteration cycle is continued 
until the error is within acceptable tolerances. After the n^* 1 step, 


0(D) 


n+1 


0(D) 


C 5T 
n n 


n 


(A2) 


where the error is given by 

&T » AT - AT 

n p p 


(A3) 


and the coefficient by 


TABLE V. - ORBIT ELEMENT AND PHYSICAL DATA FOR EARTH AND MARS 


Planet 

Perihelion 

radius, 

V 

astronomical 

units 

Eccen- 

tricity, 

e 

Perihelion 
date (ref), 
D , 

pr 

Julian 

day 

Orbital 

period, 

days 

Circular 
velocity 
at 1.1 
planet 
radii, 

Relative 
displace- 
ment of 
perihelion, 

rad 

Earth 

0.9835 

0.0167 

2444242 

365.2 

4.670 

0 

Mars 

1.3814 

.0934 

2443951 

687.5 

2.121 

4.70 


a Measured counterclockwise from Earth perihelion. 
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♦ 


n > 1 (A4) 


Since the planet orbits are nearly circular, the planet mean motions 
((M) A 2^/t) may be used to obtain good starting values which lead to very rapid 
convergence : 

= <M) AT 
0 P 

C 0 = <M) 

With the true anomaly determined in this manner, the planet radius, velocity, 
and path angle at the encounter date may be computed directly from equations (6) 
and (7). Finally, the polar travel angle p may be calculated from the true 
anomalies of Earth and Mars. For an outward (Earth-Mars) transfer, {3 ^ is the 

angle from Earth at the departure date measured counterclockwise to Mars at the 
arrival date; that is, 

Pout * a p + - W + 2jtW (A6) 

where exp is the angle from Earth's perihelion to Mars' perihelion. For an in- 
ward transfer, the angle is measured from Mars to Earth, thus 

p. = 2 jtN - a + 0-(D.) - 0 (D ) (A7) 

unw p cn cr ® v ©' v ' 

In either case, W( = 0, ±1, ±2, . . . ) is chosen to obtain the smallest positive 
value of p . 


Circular Planet Orbits 

If the planet orbits were circular, the perihelion dates and positions and 
the true anomalies used above are not well defined. In that case, it is con- 
venient to define Op = 0 and measure the planet central angles and encounter 

date relative to opposition. The iteration cycle just described then yields an 
exact result in one step; equations (A6) and (A7) reduce to 

font * % - < D d- - + ^ 

for an outward transfer, and 

%nv, - “W + ( D d- - + < A9 > 
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for inward transfer, upon noting that D in equation (Al) is replaced by 
D opp* 

These results may be used as "inputs" for the transfer trajectory calcu- 
lation shown in the Analysis section. 
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APPENDIX B 


SIMILARITY OF EARTH-MARS TRAJECTORIES IN DIFFERENT SYNODIC PERIODS 


Similarity Conditions for Earth-Mars Transfers 

For circular-coplanar planet orbits, each outward transfer possesses an in- 
ward counterpart having equal terminal radii, travel angle, and travel time when 
the Mars encounter dates are related by 

%,inw a 2D opp " %,out ( B1 ) 

as may be seen by equating (A8) and (A9) with 

3 , =* 3 • 

K out 'inl- 


and 


AT ®<f = AT^ 


The counterpart trajectories thus obtained are geometrical mirror images of one 
another and are dynamically similar. 

This mirror-image property between inward and outward transfers with equal 
travel times does not hold, in general, for the elliptic-coplanar orbits consid- 
ered herein. The property does hold, however, in one special case, when the 
lines of apsides of the two planet orbits coincide with the line of opposition. 

In this case, because of symmetry, the average planet motions are the same on 
either side of opposition, so that equations (A8) and (A9) when combined as 
above will still yield equation (Bl) . In this case, the counterpart trajectories 
are still geometrical mirror images obtained by reflection about the common line 
of opposition and apsides. They have equal terminal radii, travel angle and 
travel time, and are dynamically similar in that the AV at a given terminal 
planet is the same. 

This condition is well approximated in 1980 since Mars is close to its 
aphelion at the opposition date (2444295) and the Earth’s orbital eccentricity 
is very small compared to Mars 1 . This observation justifies the dual interpre- 
tation of the trajectory data in tables I and II as stated in the RESULTS AND 
DISCUSSION. It should be noted in addition that the dual or mirror-image prop- 
erty holds in all synodic periods, not just in 1980, as long as the reference 
date 2444295 is used. That is, an outward transfer which reaches Mars 100 days 
before the 1978 opposition would have a mirror image counterpart in an inward 
transfer departing Mars 100 days after the 1982 opposition; the two transfers 
would have equal travel times and nearly the same travel angle and terminal 
radii. 
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Effect of Synodic Period on Earth-Mars Trips 

Mars' heliocentric radius at opposition increases from 1.38 astronomical 
units in 1971 to 1.68 astronomical units in 1980, and then decreases again from 
1980 to 1989. For a given travel time and angle, the transfer requiring the 
greater radial travel will require the higher AV simply because a greater 

distance must be covered in an equal time interval. 

Round trips with N = 0 may be conveniently classified according to the 
distribution of the total available travel angle (eq. (3)) between the outward 
and return legs. If the outward travel angle 3 ou ^. is larger than P inw , the 

trip is referred to as a long-short round trip. Similarly, short-long trips 
have a long-angle return leg. 

Figure 17 shows that the choice of trajectory profile, that is, short-long 
or long-short, has a sizable effect on round trip £ AV in most synodic pe- 
riods. For either single- or double-conic trajectories, the long-short profile 
is advantageous between 1971 and 1980, while the short-long profile yields lower 
6 

AV from 1980 through 1986. This effect may be understood by considering 
the following arguments . 

CD As previously mentioned, the radius of Mars at opposition increases 
from 1.38 astronomical units in 1971 to 1.66 astronomical units in 1980, and 
then decreases again from 1980 to 1986. 

(2) For fast round trips, the angular matching requirements discussed in 
the Analysis section require the spacecraft to be at Mars within a period of 
roughly ±150 days from opposition. 

(3) In a given synodic period, the long-short profile occupies the first 
part of this ±150-day band, corresponding to early arrival at Mars. Conversely, 
the short-long profile represents late arrival in the latter part of the band. 

(4) A short-long trip is therefore likely to reach the vicinity of Mars as 
much as 300 days later than a comparable long -short trip. 

(5) The radius of Mars' orbit varies an appreciable amount during this pe- 
riod except in oppositions such as 1971 or 1980 when Mars is near one of its 
apses . 

(6) Thus, during the cycle of opposition (1971 - 1980) when the radius of 
Mars' orbit is increasing, a short-long trip, by arriving later, is required to 
travel a greater radial distance than the long-short profile. Since the angular 
velocity of Mars is decreasing during this period, the short-long trajectory 
must also describe longer central angles to match the Earth's motion. The 
short-long trip is therefore more difficult than the long-short between 1971 and 
1980. The reverse is true during the decreasing cycle between 1980 and 1989. 

For an opposition such as 1980 when Mars is near an apse, the flight profiles 
are dynamically similar and therefore require about the same ^ AV. 
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